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Abstract. Solar models through evolutionary phases of 
gravitational contraction, pre-main sequence and MS 
phases, up to current age 4.5 x 10® j/r. and 4.57 x 10® yr., 
were studied adopting different prescriptions for the equa¬ 
tion of state (EOS) and different opacity tables. 

Old EOS in Ezer’s evolutionary code that we adopted 
in previous studies includes pressure ionization, and de¬ 
generacy of electrons. In the present study the EOS is 
obtained by minimization of the free energy (Mihalas et 
al. 1990, MHD), and the radiative opacity is derived from 
recent OPAL tables, implemented by the low temperature 
tables of Alexander & Ferguson. The results are compared 
with solar models we computed with different radiative 
opacities (Cox & Stewart 1970) and different EOS, as with 
models computed by other authors. 

Finally we provide the internal run of the thermody¬ 
namic quantities of our preferred solar model which pos¬ 
sesses the following characteristics: age 4.50 x 10® yr., ini¬ 
tial He abundance by mass 0.285, parameter of the mix¬ 
ing length a = 1.82, radius and temperature at the bot¬ 
tom of the convective envelope are = 0.724 Rq and 
Tf, = 2.14 X 10® K, respectively. 


Key words: The Sun: interior-solar evolution- equation 
of state 


1. Introduction 

Recent improvements concerning radiative opacities and 
equation of state removed significant discrepancies be¬ 
tween observations and theory of stellar evolution. Up¬ 
dated solar models (Turck-Chieze et al. 1993; Charbonnel 
& Lebreton 1993; Gabriel 1994a; Basu & Thompson 1996) 
are in good agreement with the solar helioseismic data. 
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In this work models are constructed for the Sun start¬ 
ing when model becomes stable against gravitational con¬ 
traction, with different EOS and opacities, by using Ezer’s 
stellar evolutionary code (EC) (Ezer & Cameron 1967; 
Yildiz & Kiziloglu 1995). EOS of MHD (Mihalas et al. 
1990; and references therein) which uses an approach 
known as chemical picture is based on minimization of the 
free energy and is applicable for stellar modeling. We have 
chosen opacity tables of OPAL (Iglesias et al. 1992). Since 
those tables do not extend to low temperatures, Alexan¬ 
der & Ferguson opacity (1994) is used for outer layers of 
the stars. The models computed for comparison are ob¬ 
tained using Cox & Stewart (1970) (CS) opacity and EC 
EOS which takes into account degenerate electrons and 
pressure ionization in an artificial way. Since OPAL opac¬ 
ity, near the bottom of the convective zone, is larger than 
the CS opacity by a factor of about two, it strongly af¬ 
fects the structure of the convective zone. In performing 
the solar calibration, the lower luminosity due to enhance¬ 
ment of the opacity also in radiative regions is balanced 
by decreasing H (or increasing He) abundance. Thus, the 
abundance of H by mass may also be significantly altered. 

It was shown that the MHD EOS removes significantly 
the discrepancies between the observational and the the¬ 
oretical frequencies of the solar oscillations (Christensen- 
Dalsgaard et al. 1988). However we do not mention about 
the solar oscillations. The incorporation of the Coulomb 
interaction in EOS reduces both of the pressure and the 
energy. To compensate this reduction in pressure, the 
mean molecular weight per free particle, that is He abun¬ 
dance, becomes smaller, in contrast to the effect of OPAL 
opacity. The application of MHD EOS requires the calcu¬ 
lation of ionization and internal energy of each chemical 
species. Avoiding the time consuming calculation of the 
Saha equation for the heavy elements, we use two differ¬ 
ent methods, namely, Henyey (Gabriel 1994b) method and 
Gabriel & Yildiz method (1995). 

In Sect. 2 we present the basic features of the code. The 
EOS with its computational method and opacity used in 
the construction of the solar models are given in Sect 3. 
Sect. 4 is devoted to the influence of MHD EOS includ- 
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ing the results of the two methods for the ionization and 
internal energy of heavy elements. In Sect. 5 we give the 
results of the solar models with different EOS and opacity. 

2. Basics of the Modeling 

We solve the four non-linear differential equations of the 
stellar structure by the method of Newton-Raphson, from 
center all the way to the surface, with the surface bound¬ 
ary conditions given as 

L = dTTcr T^ff , +-a T^ff (1) 

where Pph is photospheric pressure, r is the optical depth 
taken to be 2/3, and the remaining symbols have the usual 
meaning. 

Mass fraction of heavy elements is taken to be 0.019. 
Chemical abundances of the solar surface found by An¬ 
ders & Grevesse (1989) are chosen. Nuclear reactions are 
from Caughlan & Fowler (1988). For convection, clas¬ 
sical mixing length theory of Bohm-Vitense (1958) and 
Schwarzschild criteria are used. 

Stellar models are employed mostly in comparative 
works. Therefore the accuracy of the code is very im¬ 
portant. For this reason the code is modified to obtain 
more accurate models. In order to acquire a higher accu¬ 
racy, two points derivative of thermodynamic quantities 
are replaced with quadratic ones. The relative changes of 
the physical variables become smaller than 10““^ for each 
shell boundaries in the solution of the equation of stellar 
structure by the Henyey method. 

3. Input Physics 

3.1. Opacity 

In order to test the influence of the uncertainities in the 
radiative opacities we have used different sets of the opac¬ 
ity. The OPAL opacity is available for the temperatures 
larger than 6000 and implemented by the Alexander 
& Ferguson opacity table at low temperatures (OA). For 
the Sun, the transition point (T = 8000 K) is chosen such 
that it always belongs to convective region in all evolu¬ 
tionary phases of the Sun. So, small discrepancies may 
not cause any problem in the solution of finite-difference 
equations. We calculate the opacity by quadratic interpo¬ 
lation in chemical composition, T and R {R = p/T^). 

Model 1 and 2 are obtained using the radiative opacity 
of CS. The opacity tables include three H-He mixtures 
in which H/He =4; 1; 0. Required opacity for the solar 
material is found by logarithmic interpolation between the 
three sets of the opacity tables. 

3.2. Equation of State: Minimization of Free Energy 

The free energy is minimized in order to obtain the Saha 
equation which is required to find the number of particles 


in each degree of ionization. This is the chemical picture 
in which atoms and molecules are described as separate 
entities. 

Since Fermi and Bose gases obey different statistics, 
their free energies, F, are different. But, disregarding 
signs, they have the same form: 

F = n + pN, (2) 

n=-PI/ = Tfcr^ln(l±e('"“"'^/''^), (3) 

i 

where summation is over all quantum states i, p, is the 
chemical potential, is the energy of this state, and k 
is the Boltzman constant. For Fermi (Bose) gas upper 
(lower) sign in Eq. (^ applies. Total number of particles, 
N, is given by the summation of the average occupation 
number over all quantum states, rij. 

3.2.1. Degeneracy of Electrons 

The free energy of a gas of particles with half integral spin 
is not as simple as that of the photon gas in those regions 
of the p — T plane where the gas is partially degenerate 
and partially relativistic. In the region of interest, there is 
no analytic solutions of the EOS except at some points in 
which both degeneracy parameter A (A = p/kT) and rel¬ 
ativistic parameter (3 {jS = kT/mc^) have extreme values. 
Therefore numerical methods should be employed for a 
semi-relativistic and semi-degenerate gas. As Qn functions 
introduced by Guess (1966) are better suited for numer¬ 
ical calculations than Fermi-Dirac functions, we express 
fie as 

TT Tnf 

^3 V (Q 4 — 4 Q 2 + 3Qo) ■ (4) 

Like fie, the number density of electrons Ue is an implicit 
function of A and (3. In terms of the functions Qn it is 
given as 

?T-e = ^ = 27r (Oa-Qi), (5) 

where A^e is the number of electrons in volume F. If A >> 
0 then perfect gas relation is recovered and A for the non- 
relativistic region can be given in terms of and T: 

e-^ = _ — _(6) 

2{2TnnkYl^ ^ ’ 

which corresponds to slightly degenerate or non¬ 
degenerate case. 

Then the energy and the pressure are given by the 
derivatives of the free energy with respect to the temper¬ 
ature and to the volume, respectively: 

= ^ (Q 4 + Q 0 ), (7) 

4 5 
TT m e 

Pe = ^ {Qi - 4 Q 2 + 3Qo). (8) 
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3.2.2. Boltzmann Gas and Partition Function 


Free energy for the Boltzmann gas can be written as fol¬ 
lows 




-fcT^A^ln 

'i',3 



f 2TrmikT\^^‘^ 

V /tO * 


.(9) 


where and Bf are number of z-type j-times ionized 
ions in volume V and partition function of these ions, 
respectively. 

Now, energy and pressure can be found for a Boltzman 
gas by differentiating Ft, given above. Since partition func¬ 
tion is also a function of temperature and density (due to 
interaction between particles), its derivatives must also be 
taken into account: 


under the influence of the Coulomb potential due to other 
charged particles in the plasma. The Coulomb interaction 
implies that slight changes occur in the energy and the 
pressure. An approximate expression for the Coulomb en¬ 
ergy is given by Landau & Lifshitz (1969) as 



1 


(13) 


where qe is proton charge and Rd is Debye radius. Ex¬ 
pression for pressure due to interaction is 


Pr = - 


3P V 


kTV 


1/2 


3/2 




mf 


(14) 


E, = kTY^Nl 



z 


E 


jm 

kT 


ainc/ 

d\nT 


As seen in Eqs. (|^ and (|^ both Ec and Pc are negative, 

( 10 ) 

since the Coulomb interaction causes a decrease in the 
energy and the pressure. 


and 


Pb 






(. d\nU^\ 

d\np j 


( 11 ) 


where n{ = /V and Bl = The second quan¬ 

tity in the parenthesis in Eq. ( 0 ) is ei^/kT where is 
the mth ionization potential of z-type ion 

Influence of interaction between atoms (ions) may 
obliterate some outer states of an atom. Electrons oc¬ 
cupying these states are not bounded to the atom any 
more. This is known as pressure ionization. MHD have 
shown how to include the effect of pressure ionization into 
the partition function based on Unsold’s theory (Clayton 
1968). According to Unsold’s theory partition function can 
be written as 


U‘ = Y,: 

k 


.-xiJkT 


( 12 ) 


where wj ^ is the survival probability of the state k for a 
j-times ionized z-type atom with j less electron, ^ and 
xl J, are the statistical weight and excitation energy of this 
state, respectively. 

If the survival probability for any state is zero, neither 
this state nor higher states contribute to the partition 
function. Apart from thermodynamical consistency the 
survival probability given by MHD has important advan¬ 
tages. The partition function doesn’t explode, and transi¬ 
tion from bound to free states is continuous since wj ^ is 
continuous. 


3.2.4. Saha Equation 


To find equilibrium number of particles it is necessary to 
minimize the free energy using stoichiometric relations de¬ 
scribing ionization (or dissociation) processes. Then the 
Saha equation is obtained with ionization potential /(, of 
Z-times ionized fc-type atom: 


K 


where 


U[ 


/ + 1 


Ui 


expjA - {li - 2{l + l)c2(iVZ2^^)'/2}/A:T 


E^ 

^,3 \ 


f din Uf d In Uj \ 

^ J 


dN\+'^ 


(15) 




(16) 


with N being the total number of particles in the system 
including electrons. 

Two remarks are called for Eq. (15) which includes de¬ 
generacy of electrons and the Coulomb interaction. Eirst, 
any additional term to the free energy appears with its 
partial derivative with respect to the number density of 
particles and as an argument of the exponential in the 
Saha equation. Second, both of the terms coming from the 
free energies of the degenerate electrons and the Coulomb 
interaction behave like processes lowering ionization po¬ 
tential. Therefore, some of the early works on EOS (Har¬ 
ris 1964; Graboske et al. 1969) treated ionization potential 
somehow as a function of the thermodynamical variables. 


3.2.3. Coulomb Interaction 

At relatively high densities atoms are very close to each 
other. A bound electron may not preserve its original state 


3.3. Computational Method of EOS 

The total energy and pressure of the plasma are the sum¬ 
mation of the partial contributions of the effects discussed 
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in the previous sections and the radiation term. Most of 
the physical quantities appearing in the equations for en¬ 
ergy and pressure of different kinds of gases can be iden¬ 
tified, explicitly or implicitly, in terms of number density 
of electrons. Therefore the EOS is calculated by iterating 
over the number of moles of electrons per unit mass, Nf. 
(Gabriel 1994b). Then the problem reduces to solve the 
Saha equation for any given value of Ng. 

We implemented the whole procedure in our stellar 
evolution code in the following way. Starting from the 
center, an initial guess of Ng is obtained by assuming the 
complete ionization, that is each species has charge Ziqg, 
which is sufficiently accurate for the center of a star. And 
then, in outer shells, the first guess of Ng is taken as its 
last computed value in the previous shell, until the surface 
is reached. 

From the value of Ng, the degeneracy parameter 
A = e~^ is found from Eq. (^. In the transition region 
0.0005 < A < 600.0, gas is partially degenerate and A 
can be recalculated from the fitting formula of Henyey 
(Gabriel 1994b; Yildiz 1996). 

Because of the low abundance of heavy elements in 
the solar mixture, probably it is unnecessary for those el¬ 
ements, to calculate the quantities such as partition func¬ 
tions and the number of elements at any degree of ioniza¬ 
tion, with high accuracy. Only for H and He very precise 
calculations need to be done. 

In order to save computer time, all that is needed is to 
calculate the effective degree of ionization (i.e., effective 
charge of each element) 




N, 


(17) 


Then, there is no need to solve the Saha equation. Two 
methods are used to compute the effective charges of ^^G, 
i4n, i6q^ |.]^g binding energies of these atoms 

(ions). The first one is Henyey’s fitting method (Gabriel 
1994b; Yildiz 1996) which uses NgkT and the chemical 
potential of electron ^g as input parameters. The second 
method is due to Gabriel & Yildiz (1995), which is similar 
to Henyey’s method but more precise. It takes into account 
Debye radius as an additional parameter. 

Since heavy elements are rare in the atmospheres of 
the stars only the molecules of hydrogen, H 2 , HJ, and H“ 
are considered. The abundance and the energies of these 
molecules are solved by the method of Vardya (1961). 


4. Influence of MHD EOS 

In order to compare the results and the influence of MHD 
EOS, the required density, temperature and chemical com¬ 
position are taken from a model of the present Sun (Model 
A, see Table 1). Henyey method is used for ionization of 
heavy elements, unless stated otherwise. 

Contribution of the Coulomb interaction, the degener¬ 
acy of the electrons and derivative of the partition function 


Fig. 1. The fractional contribution to the pressure due to 
Coulomb interaction, degeneracy of electrons and derivative 
of the partition functions of H, He, and He^ with respect to 
density, is plotted as a function of the logarithm of the tem¬ 
perature for a selected solar model (case A below). 


to the total pressure can be seen in Fig.l. For the Model 
A P/{P — P') is plotted with respect to log(T) of tem¬ 
perature from surface to the center throughout the solar 
model, where P is the total pressure, and P' is any of frac¬ 
tional contribution to the pressure. Thin solid line is for 
the Coulomb pressure. The ratio is less than or equal to 
one, since the Coulomb term causes a decrease in pressure. 
At the surface the Coulomb pressure is zero. As the tem¬ 
perature and the density, and therefore the degree of ion¬ 
ization increase, fractional change becomes a few percent. 
Contribution of the Coulomb potential to the total pres¬ 
sure is maximum (6%) between log(r) = 4.6 and 4.7, 
at which point the most abundant species (hydrogen and 
helium) are completely and singly ionized, respectively. 
This means that the Debye radius is sufficiently small. As 
the temperature and the density increase the contribution 
becomes smaller, until partial degeneracy of the electrons 
sets in near the center of the Sun. 

The pressure of the degenerate-electrons calculated 
by Henyey fitting formula (also calculated by method of 
Yildiz & Eryurt-Ezer (1992) using Guess functions, the 
difference is less than 0.01%), represented by dots, de¬ 
viates from one at the central part. Here NgkT is sub¬ 
tracted from the total electron pressure. At the center of 
the Sun degeneracy parameter (A) is equal to —1.559 
which implies a slight degeneracy, and the contribution of 
degenerate-electrons is about 2%. 

Thick solid line in Fig.l. is a superposition of three 
Gaussian curves. From left to right, the first one is dom¬ 
inantly the density derivative of hydrogen partition func¬ 
tion which is maximum (w 0.5%) at log(r) w 4.5. Sec- 
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Fig. 2. The run of the effective charge of several elements as 
a function of the logarithm of the temperature in the model of 
Fig. 1. The Saha equation is solved for H and He, and Henyey’s 
method is used for heavy elements. 


Fig. 3. The run of the effective charge of several elements as 
a function of the logarithm of the temperature in the model 
of Fig. 1. The Saha equation is solved for H and He, and the 
method of Gabriel & Yildiz (1995) is used for heavy elements. 

ond one is due to the density derivative of He partition 
function, which is maximum at about log(T) « 4.75 and 
its contribution to the total pressure is less than 1%. Last 
one is due to the density derivative of He’*' partition func¬ 
tion, which is non-zero in an important part of the solar 
interior but with a lower contribution to the total pressure 
than the others. The base temperature of the convective 
zone in the Model A is about 2.14 x 10® K. This corre¬ 
sponds to a point (log(T) = 6.3) at which He+ starts to 
appear. 


Fig. 4. The ratio of adiabatic gradients resulting from MHD 
EOS with Gabriel & Yildiz (1995) (solid line) and with Henyey 
(dots) methods for ionization of heavier elements, to that ob¬ 
tained adopting the EOS of old routine in Ezer’s code (EC) for 
the selected model of Fig. 1. 

The Saha equation, Eq (15), is solved exactly for H 
and He. Results obtained from Henyey method, the effec¬ 
tive charges divided by the atomic numbers of each ele¬ 
ment, throughout the Model A, are given in Fig.2. Due to 
its low ionization potential hydrogen (solid line with as¬ 
terisks) is rapidly completely ionized (at log(T) « 5.2). 
For the effective charge of He (solid line with circle) two 
phases are seen. Even at the center, none of the heav¬ 
ier elements, C (solid line), N (dashed line), O (solid line 
with triangle), Ne (dashed line with asterisks) are com¬ 
pletely ionized. Owing to rough calculation, each line of 
the heavy elements cross each other near the center. On 
the other hand, at the center the radiation and the gas 
pressure (and energy) are sufficiently large, therefore the 
computed values are within an acceptable range. This is 
not the case at the surface where binding energy has an 
important role in the total energy. While H and He are 
neutral, each of the heavier elements is in the state of first 
ionization. 

A physically more reliable method is that of Gabriel 
& Yildiz (1995). As seen in Fig.3, at the middle of the 
Sun, ionization degrees of the heavier elements are a little 
bit smaller than, but comparable to those given in Fig.2. 
At the center and at the surface, the values of ionization 
degree are not in contradiction with the basic truths of 
atomic physics. All atoms are neutral at the surface. Ion¬ 
ization of C starts before H, and ionization of N, 0,and Ne 
starts after H but before He, as expected from the atomic 
physics. 

The effect of the MHD EOS on the adiabatic gradient 
is shown in Fig.4. There we plot the ratio of the adiabatic 
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Table 1. Comparison of the present solar models with different EOS, opacities and data of the Sun. 


Models 

Y 

a 

Rc/Rq 

Ti/W° 

Teff 

Tc/lO” 

Pc 

W 

Pc/10"' 

1 (CS EC Dl) 

0.224 

1.32 

0.19 

1.24 

5804 

15.00 

136.4 

0.438 

2.244 

2 (CS MHD-tH Dl) 

0.205 

1.46 

0.20 

1.31 

5794 

14.91 

132.0 

0.478 

2.147 

3 (OA EC Dl) 

0.302 

1.57 

0.27 

2.07 

5799 

15.80 

151.9 

0.324 

2.342 

4 (OA MHD+H Dl) 

0.284 

1.76 

0.27 

2.11 

5797 

15.70 

149.3 

0.338 

2.271 

A (OA MHD+GY Dl) 

0.285 

1.82 

0.276 

2.14 

5799 

15.70 

147.1 

0.345 

2.269 

B (OA MHD+H Dl) 

0.285 

1.78 

0.277 

2.15 

5799 

15.69 

147.5 

0.343 

2.267 

C (OA MHD+GY D2) 

0.283 

1.82 

0.276 

2.16 

5792 

15.67 

148.6 

0.339 

2.270 

GL (1993) 

0.278 

1.66 

0.281 

2.13 


15.56 

150.7 



BU (1988) 

0.271 


0.277 

2.11 


15.62 

148 

0.341 

2.290 

TGGD (1988) 

0.275 

2.11 

0.270 

2.04 


15.52 

147.2 



GGK (1989) 

0.291 

1.89 

0.286 

2.29 

5770 

15.68 

162.4 

0.352 

2.2 78 

SBF (1990) 

0.278 

2.07 

0.260 

1.96 


15.42 

146.6 



LD (1988) 

0.278 

2.16 

0.265 

2.02 

5781 

15.54 

148.0 

0.352 

2.278 

GDR (1996) 

0.265 

1.77 

0.274 

2.10 


15.45 

147.2 

0.372 



gradient resulting from the Gabriel & Yildiz (solid line) 
and from the Henyey method (dotted line) to the value 
obtained with EC EOS are given. Both are larger than 
the old gradient, and three peaks correspond ionization 
zones of H, He and He"*". 

5. Results and Discussion 

Recently we have gathered very important information 
about the internal structure of the Sun by helioseismolog- 
ical investigations. Distance of the convective zone from 
the center of the Sun is determined from p-mode oscilla¬ 
tions as 0.713 ± 0.003 Rq by Christensen-Dalsgaard et al. 
(1991). Photospheric abundance of He by mass is found 
as 0.242 ± 0.003 (Hernandez & Christensen-Dalsgaard 
1994). There are several attempts to determine the cen¬ 
tral density of the Sun from helioseismological data. Re¬ 
sults obtained by Dziembowski et al. (1994) are in between 
130—150 g cm~^, while Gough & Kosovichev (1988) quote 
w 160 g cm~^ and Vorontsov & Shibahashi (1991) quote 
a range of 110 — 115 p cm~^. 

Two sets of the solar data are used in the calculations. 
In one set (Dl), the age of the Sun is taken as 4.5 x 10® yr., 
and its luminosity and radius are 3.90 x 10^® erg s~^ and 
6.951 X 10^® cm, respectively. In the other set (D2) (case 
C below), the luminosity and the age of the Sun are 
taken from Bahcall et al. (1995), namely Lq = 3.8440 x 
10®® erg s~^ and t© = 4.57 x 10® pr. The radius of the 
Sun is Rq = 6.9598 x 10®® cm (Sackmann et al. 1993). In 
both of the sets the solar mass is 1.985 x 10®® g. 

For different combinations of EOS’s and opacities and 
for different solar data, starting from threshold of stabil¬ 
ity point at which gravitational and internal energies are 
nearly the same, we constructed a series of evolutionary 
models for the Sun. Each model given in Table 1 is ob¬ 
tained by changing initial mass fraction of hydrogen, he¬ 
lium and convective parameter a, in order to fit the Sun’s 


luminosity and radius to the present data, with different 
accuracies. The first column in Table 1 represents the ini¬ 
tial abundance of He by mass, and a is given in the second 
column. The depth of the convective zone in terms of so¬ 
lar radius, RdR q, is shown in the third column, and base 
temperature of the convective zone, T^, is in the fourth 
column. Last five columns are the effective temperature 
(Teff), the central temperature (Tc), density(pc )5 abun¬ 
dance of hydrogen by mass at the center (Xc), and the 
central pressure (Pc), respectively (in cgs). All the mod¬ 
els, except Model C, is with Dl. The first three models 
are constructed with an accuracy of a few percent, and the 
fourth model with an accuracy of 0.5%. Models A, B and C 
are more precise models (accuracy is 10“^). The difference 
between the Model A and the Model B is the method of 
heavy element ionization. The latter (also Model 4) uses 
Henyey’s method (MHD-I-H), the former is obtained by 
the method of Gabriel & Yildiz (MHD-fGY) as the Model 
C. 

When we compare each of the first four models with 
the model of the same EOS but different opacity (Model 
1 with 3, and Model 2 with 4), it is seen that OA opacities 
strongly influence the structure of the Sun. The enhance¬ 
ment of the opacity (OPAL) below the convective zone en¬ 
larges the zone’s size toward the center of the Sun. When 
MHD is used, OA opacities enlarge the convective zone by 
about 40%, and by 35% if EC EOS is used. Increase in base 
temperature of convective zone is about 61% with MHD 
EOS and 67% with EC EOS whenever OA opacities are 
employed in place of CS opacity. While He abundance by 
mass is exceeding the helioseismological result, the central 
density is about 150 g cm~^, and the base temperature of 
the convective zone is larger than 2 x 10® K in the mod¬ 
els with OA. Because of increase in opacity, less energy 
reaches the surface. This is compensated by an increase in 
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Table 2. The present Sun model (Model A) with OPAL opacity and MHD EOS by using method of Gabriel and Yildiz for 
ionization of heavy elements 


Sh.n 

M/Mq 

R/Rq 

P 

T 

P 

E 

L/Lq 

H 

1 

0.0000 

0.000 

.1471E-f3 

.1570E-h8 

.2269E-fl8 

.2289E-hl6 

0.000 

.3455 

7 

0.0004 

0.016 

.1446E-f3 

.1560E-h8 

.2232E-fl8 

.2290E-hl6 

0.004 

.3518 

17 

0.0017 

0.019 

.1436E-f3 

.1557E-h8 

.2217E-fl8 

.2291E-hl6 

0.007 

.3543 

41 

0.0022 

0.028 

.1397E-f3 

.1542E-h8 

.2158E-fl8 

.2292E-M6 

0.021 

.3640 

66 

0.0074 

0.042 

.1309E-f3 

.1510E-h8 

.2028E-fl8 

.2301E-hl6 

0.066 

.3889 

106 

0.0550 

0.089 

.9411E-f2 

.1355E-h8 

.1470E-fl8 

.2321E-hl6 

0.368 

.5126 

126 

0.2019 

0.155 

.5372E-f2 

.1100E-h8 

.7553E-fl7 

.2089E-hl6 

0.804 

.6373 

133 

0.2807 

0.182 

.4166E-f2 

.1002E-h8 

.5442E-fl7 

.1941E-hl6 

0.901 

.6636 

155 

0.5582 

0.279 

.1539E-f2 

.7289E-h7 

.1490E-fl7 

.1438E-hl6 

0.996 

.6921 

179 

0.7825 

0.393 

.4271E-fl 

.5241E-h7 

.2972E-fl6 

.1036E-hl6 

1.000 

.6959 

219 

0.9424 

0.593 

.5375E-fO 

.3195E-h7 

.2281E-fl5 

.6330E-hl5 

1.000 

.6962 

236 

0.9701 

0.683 

.2346E-fO 

.2511E-h7 

.7823E-fl4 

.4977E-hl5 

1.000 

.6962 

256 

0.9787 

0.724 

.1659E-fO 

.2144E-h7 

.4720E-fl4 

.4245E-hl5 

1.000 

.6962 

280 

0.9912 

0.808 

.8149E-1 

.1340E-h7 

.1444E-fl4 

.2634E-hl5 

1.000 

.6962 

298 

0.9958 

0.857 

.4752E-1 

.9381E-h6 

.5878E-fl3 

.1831E-hl5 

1.000 

.6962 

356 

0.9994 

0.933 

.1265E-1 

.3927E-h6 

.6457E-fl2 

.7441E-hl4 

1.000 

.6962 

369 

0.9996 

0.945 

.9133E-2 

.3174E-h6 

.3747E-fl2 

.5947E-hl4 

1.000 

.6962 

389 

0.9998 

0.961 

.5046E-2 

.2167E-h6 

.1394E-fl2 

.3937E-hl4 

1.000 

.6962 

409 

0.9999 

0.970 

.3053E-2 

.1589E-h6 

.6090E-fll 

.2753E-hl4 

1.000 

.6962 

425 

0.9999 

0.976 

.1993E-2 

.1246E-h6 

.3067E-fll 

.2005E-hl4 

1.000 

.6962 

442 

1.0000 

0.981 

.1345E-2 

.1012E-h6 

.1656E-fll 

.1463E-hl4 

1.000 

.6962 

458 

1.0000 

0.985 

.8592E-3 

.8004E-h5 

.8222E-flO 

.9747E-hl3 

1.000 

.6962 

475 

1.0000 

0.987 

.5616E-3 

.6345E-h5 

.4179E-flO 

.6144E-hl3 

1.000 

.6962 

493 

1.0000 

0.990 

.3518E-3 

.4986E-h5 

.2000E-fl0 

.3096E-hl3 

1.000 

.6962 

512 

1.0000 

0.992 

.2014E-3 

.3948E-h5 

.8721E-f09 

.3266E-hl2 

1.000 

.6962 

522 

1.0000 

0.993 

.1454E-3 

.3527E-h5 

.5502E-f09 

-.9956E-hl2 

1.000 

.6962 

533 

1.0000 

0.994 

.1043E-3 

.3179E-h5 

.3475E-f09 

-.2183E-hl3 

1.000 

.6962 

543 

1.0000 

0.995 

.7647E-4 

.2910E-h5 

.2282E-f09 

-.3169E-hl3 

1.000 

.6962 

563 

1.0000 

0.996 

.3947E-4 

.2464E-h5 

.9505E-f08 

-.4962E-hl3 

1.000 

.6962 

583 

1.0000 

0.997 

.2024E-4 

.2139E-h5 

.4031E-f08 

-.6460E-hl3 

1.000 

.6962 

594 

1.0000 

0.997 

.1433E-4 

.2002E-h5 

.2607E-f08 

-.7143E-hl3 

1.000 

.6962 

603 

1.0000 

0.998 

.1045E-4 

.1892E-h5 

.1756E-f08 

-.7723E-hl3 

1.000 

.6962 

613 

1.0000 

0.998 

.7633E-5 

.1792E-h5 

.1188E-f08 

-.8267E-hl3 

1.000 

.6962 

623 

1.0000 

0.998 

.5036E-5 

.1674E-h5 

.7103E-f07 

-.8944E-hl3 

1.000 

.6962 

643 

1.0000 

0.999 

.2610E-5 

.1508E-h5 

.3158E-f07 

-.9943E-hl3 

1.000 

.6962 

663 

1.0000 

0.999 

.1296E-5 

.1348E-h5 

.1325E-f07 

-.1096E-hl4 

1.000 

.6962 

674 

1.0000 

1.000 

.8965E-6 

.1266E-h5 

.8333E-f06 

-.1150E-hl4 

1.000 

.6962 

694 

1.0000 

1.000 

.4730E-6 

.1108E-h5 

.3607E-f06 

-.1250E-hl4 

1.000 

.6962 

705 

1.0000 

1.000 

.3356E-6 

.9903E-h4 

.2195E-f06 

-.1311E-hl4 

1.000 

.6962 

715 

1.0000 

1.000 

.2679E-6 

.8210E-h4 

.1407E-f06 

-.1363E-hl4 

1.000 

.6962 

716 

1.0000 

1.000 

.2658E-6 

.8030E-h4 

.1363E-f06 

-.1366E-hl4 

1.000 

.6962 

717 

1.0000 

1.000 

.2642E-6 

.7829E-h4 

.1320E-f06 

-.1370E-hl4 

1.000 

.6962 

727 

1.0000 

1.000 

.2421E-6 

.6432E-h4 

.9894E-f05 

-.1387E-hl4 

1.000 

.6962 

739 

1.0000 

1.000 

.1682E-6 

.5799E-h4 

.6197E-f05 

-.1393E-hl4 

1.000 

.6962 


He abundance. The convective parameter alpha increases 
by about 30%, when OA is employed in place of CS. 

If we compare Model 1 with 2, and Model 3 with 4, it 
becomes obvious that the depth and the base temperature 
of the convective zone is not very sensitive to the EOS. The 
values of He abundance and a are changed by EOS. In¬ 
corporation of the Coulomb interaction reduces the energy 
and pressure by some fraction. In order to compensate this 


reduction in pressure, mass fraction of He (or molecular 
weight per free particle) decreases by about 6% in case 
of CS, and 10% in case of OA opacities. The only differ¬ 
ence between the ways that Model A (Gabriel & Yildiz) 
and B (Henyey) are obtained, is the calculation method 
of the ionization of heavy elements and their internal en¬ 
ergies. There are small differences between these models. 
The convective parameter a of Model A is a little bit larger 
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than that of Model B, since the pressure scale height (the 
number density of electrons) of Model A in outer layers is 
less than that of Model B. 

The small differences between Models B and 4 are due 
to low accuracy of Model 4. The Models A (Dl) and C 
(D2) have different solar data. The Model C has lower 
luminosity and higher age than the Model A. 

In Table 1, we also give models constructed by Char- 
bonel & Lebreton (CL) (1993), Bahcall & Ulrich (BU) 
(1988), Turck-Chieze et al. (TCCD) (1988), Cox et al. 
(CGK) (1989) , Sackmann et al. (SBF) (1990), Lebreton 
& Dappen (LD) (1988), and Ciacio et al. (CDR)(1996). 
Similar to our models, all the models obtained by differ¬ 
ent authors are standard models, that is there is no rota¬ 
tion, no diffusion process and magnetic field is negligible. 
Except model of CGK, the initial mass fraction of He in 
Models A, B and C is a little bit larger than those found 
by other researchers. The extreme model of CGK, which 
has the largest value of He mass fraction, base tempera¬ 
ture of the convective zone, and the central temperature, 
has an age of 4.66 x 10® yr. The EOS of CGK is very sim¬ 
ilar to EC EOS, and Iben’s fitting formula (1975) is used 
for opacity. The model of CDR which has the same solar 
data as Model C except the solar mass and heavy element 
abundance, is obtained by using both the EOS and the 
opacity of OPAL. The significant difference between these 
two models is in He abundance due to their data of high 
solar mass and low metallicity. The other models take the 
age of the Sun as 4.6 x 10® yr. Therefore, their He mass 
fractions are closer to, but less than, that of our models 
with OA and MHD EOS. Almost all the models use the 
same nuclear reaction rates of Caughlan & Fowler. Model 
of CL, which uses OPAL opacity and tables of MHD EOS, 
is close to our best models. There is a small difference in 
the value of a (the difference in He mass fraction is due 
to different age of the Sun). This difference stems possi¬ 
bly from their different heavy element abundance and low 
opacity table. They use different initial mixture for heavy 
elements. Model of BU includes the Coulomb interaction 
in EOS and the Los Alamos Opacity (Cox et al. 1991). 

The physical variables which determine the structure 
of the Sun is summarized in Table 2. These values (in cgs) 
are from the Model A which has the accuracy of 10“"'’ and 
is obtained by using MHD-I-GY and recent opacities. Eor 
accuracy of the model, number of shells is increased to 739, 
during the evolution. The energy that the Sun radiates is 
produced within the inner 0.25 Mq core. Convective zone 
is between shells 256 and 727, where radiative temperature 
gradient exceeds adiabatic temperature gradient. Its mass 
is only 2% of the total mass of the Sun, and its distance to 
the center is O.724i?0. The total internal energy in outer 
shells is negative, since binding energy of an atom having 
electronic configuration is considered negative. 

Our conclusion is that solar models obtained by MHD 
EOS and OPAL opacities are in closer agreement with 
the results the helioseismology. The remaining small dif¬ 


ferences can be removed by taking into consideration the 
diffusion process of He and heavy elements. With MHD 
EOS and OPAL opacity, and with diffusion process, Basu 
& Thompson (BT) (1996) gives the distance from bottom 
of the convective zone to the center as 0.714 Rq. When 
one incorporates diffusion processes, Bahcall et al. (1995) 
emphasize that the decrease in mass fraction of He is about 
10% which gives better agreement with the observation. 
From Table 1, one sees that the main effect upon the initial 
He abundance is provided by the adoption of the recent 
opacities which tend to increase the abundance, while the 
MHD EOS has a smaller effect in the opposite direction 
Structure of the Sun is not very sensitive to the meth¬ 
ods of Henyey (Model B) and Gabriel & Yildiz (Model A) 
for the ionization of heavy elements. But the method of 
Gabriel & Yildiz is better than Henyey method for rapid 
fitting processes, and, near the surface of the Sun, it gives 
neutral heavy elements as expected. 
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